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We present experimental and theoretical results on hole-spin dynamics in bulk GaAs after ultra- 
' fast optical excitation. The experimental differential transmission are compared with a dynamical 

, calculation of the momentum-resolved hole distributions, which includes the carrier-carrier, carrier- 

phonon and carrier-impurity interaction at the level of Boltzmann scattering integrals. We obtain 
good agreement with the experimentally determined hole-spin relaxation times, but point out that 

Q depending on how the spin-polarization dynamics is extracted, deviations from an exponential de- 
cay at short times occur. We also study theoretically the behavior of the spin-relaxation for heavily 
QQ ' p-doped GaAs at low temperatures. 
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Research on spin dynamics and spintronics in semiconductors and metals covers an immense variety of applications 
in information storage and manipulation in carrier spins in solid state systems>i"— Both storage and manipulation of 
00 , spins are limited by spin relaxation, so that an understanding of spin relaxation processes is important."*'^ For ultrafast 
spin dependent dynamics on picosecond timescales, simplified relaxation-time approximations are no longer justified 
and a fully microscopic understanding of the complex spin dynamics is needed<^"— While electron spin-dynamics 
continue to be extensively investigated in semiconductors due to their extremely long spin-relaxation times,- hole- 
spin dynamics are intriguing for different reasons. In III-V semiconductors, hole dynamics is inherently different 
from that of electrons because of the strong spin-orbit coupling of the hole states. Thus, hole-spin relaxation and 
momentum (energy) relaxation occur on a comparable ultrashort timescale. 

In this paper, we connect differential transmission measurements of the hole-spin dynamics^'' with a microscopic 
' O ' calculation that includes the anisotropic band-structure as well as carrier-carrier and carrier- phonon/impurity scat- 
C . tering mechanisms at the level of Boltzmann scattering integrals. Although our approach does not take into account 
^ all the intricacies of the dynamical multi-band screening (including phonon-plasmon coupling)— and interband 
polarizations, it captures essential aspects of the spin dynamics of holesJ^ We show that this calculation compares 
well with measured hole-spin relaxation times, but that deviations from an exponential spin decay occur for short 
times depending on how the relaxation time is extracted. Further, we present results for the hole-spin relaxation at 
high p-doping, which underlies a recent treatment of GaAs-based ferromagnetic semiconductors!^ 

Recent theoretical treatments of hole spin (and charge) dynamics in bulk GaAs have focused on the calculation of 
energy-dependent spin relaxation-rates due to phonon/impurity scattering^ the investigation of coherent hole (spin) 
■ dynamics in GaAs,^'''^^ and the influence of Dyakonov-Perel type precession effects due to spin-dependent splittings 
, between hole bands in GaAs quantum wellsJ^ Very recently, additional contributions due to "small" band structure 
effects and the spin-orbit contribution to the hole-phonon interaction for holes in bulk GaAs have also been analyzedf^^ 
From an experimental point of view, hole-spin polarizations can be created efficiently using optical orientation 
techniquesJS Because of the ultrashort lifetimes^^ involved, it is difficult to unambiguously study hole dynamics 
in degenerate pump-probe experiments (i.e., experiments with the same pump and probe wavelengths) due to the 
competing presence of the coherent artifact on the same timescale^S. and, to a lesser degree, the electron spin dynamicsi^ 
Hole spin study, therefore, requires measurement of the spin-dependent carrier dynamics on ultrashort timescales 
using a time-delayed, mid- infrared probe beam (i.e., a nondegenerate pump-probe experiment), which was realized 
only much later"'^^ than similar experiments on electrons<^ 

We focus here on the dynamics of the spin-polarized heavy hole subband. We have performed nondegenerate, 
polarization-resolved pump-probe spectroscopy at 300 K to measure the hole spin relaxation lifetime in intrinsic 
GaAs as a function of input pump fluence, J-. In this setup, 1 W of average power from an 80 fs titanium:sapphire 
laser (Spectra Physics Tsunami) operating at an 80 MHz repetition rate is used to pump a custom-designed optical 
parametric oscillator based on periodically-poled lithium niobate, as described in further detail in ref.[23l A portion of 
the titanium:sapphire beam (< 300 mW) is circularly polarized using a quarter wave plate and is used to photoexcite 
spin-polarized electron-hole pairs in the GaAs sample. The pump beam induces a circular birefringence and results 
in a time-dependent rotation of the transmitted probe polarization. To isolate the dynamics of the hole spin state 
from the electron spin, we probe the occupancy of the heavy hole and/or light hole subbands using a synchronized 
mid-infrared idler pulse (A = 3200 nm), which probes the induced transitions from the split off subband into the 
heavy hole subband. This wavelength directly probes states from the lower split-off hole subband to the heavy hole 
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FIG. 1. Fluence dependence of the time-resolved differential transmission, AT/T{t), at 300 K. 



subband at the same quasimomentum (k « 0.3 nm~^) as the initial photoexcitation at 800 nm. 

Figure [1] shows the measured field rotation, AT/T(t), of the mid-IR probe from — 2.6 to 13.8 /iJ cm^^. The 
input probe beam is linearly polarized (£), which is an equal admixture of both the co-circularly polarized ((t+) and 
counter circularly polarized(i7_). A MgF2 linear polarizer is rotated 90° with respect to the input polarization (i.e., 
along y) and is used to measure the field rotation of the transmitted probe, which depends on the time-dynamics of 
the the heavy hole occupancy state near the quasimomentum, fc, of photogeneration. Using the measured absorption 
coefficient of GaAs,- - the resulting carrier concentrations range from 5 x 10^^ cm~^ = 2.6 /.tJ cm^^) to 2 x 10^^ cm~^ 
(J- — 13.8 /iJcm"^). We fit the measured field rotation to a single exponential model and account for the finite pump 
and probe pulse widths by convolving the model function with a gaussian with a width of 120 fs (see ref. |10| for more 
details of the fitting model). We find no significant trend in the extracted lifetimes, which range from 120 fs to 145 fs, 
within the experimental error (±20%) of this experiment. 

The calculation of the hole-spin dynamics follow the procedure outlined in ref. [l^ The electron and hole states 
around the fundamental bandgap needed as input for the dynamical calculation are determined from an 8-band Kane 
model H{k) with six hole and two electron bands containing terms up to second order in fc.^^ By diagonalizing 'H(fc) 
we obtain the "intelligent basis" (in the sense of ref. [26h for the hole states: ji^, k) and their energy eigenvalues 
The label ly = {b,p) includes the band index & = E, HH, LH, SOH, (for electrons, heavy holes, light holes, and split-off 
holes) and the pseudospin p = 1,2. The pseudospin can be introduced because the quasiparticle dispersions of all 
four types of carriers are (nearly) doubly degenerate. 

The dynamical equations governing the time evolution of the incoherent carrier distributions, ^, under the 
influence of optical fields and scattering are>i^^ 

r: = r"-(l-n ,-) - r°"in r (1) 

i^.k^ t^'^' v,k v,k y I 

(plus a contribution from the optical excitation, which neglects the hole band coherences between hole states of the 
"intelligent basis," or, equivalently, only describes contributions to the hole spin relaxation of the Elliott type^^ that 
arise from the k dependent mixing of different spins in the states \v^k). Contributions to the spin relaxation via 
coherent spin precession, i.e., Dyakonov-Perel-type spin dynamics;^ are neglected. 

The dynamical in-scattering rate consists of the carrier-carrier interaction contribution (i.e., eq. (2) in Ref. [11) and 
the carrier-phonon interaction 

^1"^-? = Y E |M,,An(^i,^+<7>,fc)Pv,fe+,-[(l + ^?.A)<5(Ai?_) + 7V,,A<5(Ai?+)] . (2) 

where Ai?± = Sy,k — s^^ \k+^ '^^q,\- There is also a contribution similar to eq. ^ due to carrier-impurity scatteringii 
Here, -^^9, a 7 and iOq_\ are the phonon occupation numbers, carrier-phonon coupling matrix elements, and phonon 

dispersions for LO and LA phonons, respectively.^^ The out-scattering rates r°"* are obtained from F'" by exchanging 
(1 — n) with n and (1 + A^) with N . Equation ^ and its counterpart for Coulomb scattering describes two-particle 
scattering processes connecting states \v^k) — > \v\^k ^ q) (and \v2^k\ ^ q) — > \vz^k\)) with different average spin. 
The pronounced anisotropy of the single-particle states, which is due to the spin-orbit interaction, is included in 
our calculation in the overlaps and the carrier energies, which enter the energy-conserving delta functions in the 
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FIG. 2. Theoretical (curve) and experimental values (symbols) for fluence dependence of hole-spin relaxation in intrinsic GaAs 
at room temperature. 



scattering rates. The anisotropies and high-dimensional scattering integrals result in a challenging numerical problem. 
We incorporate the effects of the anisotropy by expanding into spherical harmonics Yi^m{k") up to order ^ = 4 

and retain only the expansion coefficients with radial or cubic symmetry, as these are the dominant symmetries of 
the Hamiltonian "H.^^ This procedure, together with the use of static screening, w'*(g) = e^/[eoebg('?^ + k^)], reduces 
the numerical complexity. Additional simplifications are achieved by including only the heavy-hole bands in the 
dynamical calculation as the number of optically excited light holes is much smaller, and by replacing the non- 
equilibrium electron distributions by equilibrium distributions with the lattice temperature. We then calculate the 
spin polarization of the HH band. In ref. [l^, it was shown that the spin polarization is not identical to the DT signal 
that is measured, but the deviations are small enough for the accuracy of the present comparison. 

We include the carrier excitation by an ultrashort optical pulse using a (T+ polarized plane wave traveling in the 
(001) direction (z) to generate an initial condition for the carrier distributions 

= 0) = E I • E\h{rv^ - - • (3) 

Here, dp,^{k) = e(/i, fclrji^, fc) are the dipole-matrix elements, E is the laser field, and hio is the photon energy of the 
exciting field. A Gaussian broadening function g{e^ peaked at e — (^hu ~ £p £ ~ £^ accounts for the spectral width 
of the pulse (15meV). 

To compare with experimental results in Fig. [21 we convert the experimental laser fluence into excited carrier 
densities, and extract the relaxation time from an exponential fit to the spin polarization at the probe-laser wavelength. 
The experimental results are in good agreement with our calculations and a recent theoretical study. Fig. [2] shows 
that hole spin-relaxation is rather insensitive to excited carrier densities on the order of 10^^ cm^'^. The calculation 
also predicts a very weak temperature dependence with a hole spin relaxation time of 250 fs for an excited density of 
lO^^'cm"^ at 4K, which is somewhat shorter than the result predicted by the Elliott- Yafet spin relaxation-time.^^ 

In Fig. |3l we simulate the momentum dependence of the spin polarization dynamics. We model the excitation by 
an 800 nm pump pulse that excites a carrier density of lO^^cm"^. This leads to anisotropic initial hole distributions 
that are peaked at momentum k = 0.3 nm~^ (after integration over the angular variables). Fig.[3]shows different spin 
polarizations calculated using the angle-averaged distributions n{k,t) for different modulus of the hole momentum k. 
The decay of the spin polarization near the peak position of the initial distribution is nearly exponential, as can be 
seen from the linear curve shape on the logarithmic scale. The evolution of the spin polarization for larger and smaller 
momenta (i.e, away from the momenta of the initially photoexcited holes) is significantly faster and non-exponential 
for the initial ^ 100 fs. After the initial non-exponential dynamics the spin-relaxation becomes exponential again with 
the same time constant. The spin relaxation-time therefore depends on the way it is extracted from the calculation, 
as opposed to the case of electron spin-relaxation in p-doped GaAs where the spin relaxation time can be quite 
rigorously defined. The dashed lines in Fig. [3] correspond to experimental situations where pump and probe pulses 
are detuned, as was investigated in ref. [13 (see Fig. 5 therein) , where no change in spin relaxation times for different 
detuning between pump and probe pulses was found. The apparent discrepancy likely results from the analysis of the 
experimental data via a convolution of Gaussians, which are used to describe both the pulse shapes and the dynamics 
of the material response. The analysis of the pump-probe experiment therefore presupposed an exponential decay of 
the polarization, which turns out to be only somewhat longer than the experimental time resolution of about lOOfs^Si. 
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FIG. 3. Computed hole-spin polarization dynamics at different hole momenta for fixed excitation at k = 0.3 nm ^ in intrinsic 
bulk GaAs. The excited carrier density is n = 10^ cm~^. 
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FIG. 4. Computed p-doping dependence of hole-spin relaxation in GaAs at low temperatures (solid line). The dashed line 
results if only scattering by from ionized donors is included. 

For degenerate pump and probe, this analysis yields the exponential decay of the spin polarization in agreement with 
the calculated result in Fig. [3] For non-degenerate situations the non-exponential features in Fig. [3] on timescales 
below 100 fs are not resolved. Their resolution would require a time resolution well below 100 fs. 

Finally, we examine the hole-spin relaxation dynamics in heavily p-doped GaAs at low temperatures (5 K) theo- 
retically. By choosing these parameters, we estimate the hole-spin relaxation time for GaMnAs in the ferromagnetic 
phase at low temperatures, which has similar carrier doping levels. We model the optical excitation of nonequilibrium 
carriers in p-doped GaAs by a pump laser with a wavelength of 400 nm. Choosing this wavelength ensures that the 
optical transitions are not blocked by holes even for very high doping densities. We assume initial carrier distributions 
of the form n° = 71^°^ + n'jf'^, which contain both the influence of the doping and the optical excitation. All itinerant 
holes introduced by the p-doping are assumed to be thermalized at the lattice temperature; they are modeled by 
unpolarized Fermi-Dirac electron distributions n'^°P = /(e|,) at lattice temperature with carrier density equal to the 
density of dopants. For the optical excitation we assume again an ultrashort pulse with a width of 15 meV. Fig. 4 
shows that the hole-spin relaxation drops to approximately 10 fs for doping densities typical for GaAs based magnetic 
semiconductors. At these high doping densities, the rapid momentum relaxation due to scattering from the ionized 
donors provides the main relaxation mechanism. The assumption of a hole spin relaxation time of '^10 fs in ref. H 
makes it possible to neglect a spin-bottleneck effect for ferromagnetic GaMnAs. Our results, based on a dynamical 
calculation, support this assumption. 
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In conclusion we examined ultrafast hole-spin relaxation using differential transmission measurements and a mi- 
croscopic theoretical approach for different excitation conditions. Theory and experiment show that the hole-spin 
relaxation time is approximately 100 fs and independent of moderate changes of the pump fluence at room temper- 
ature. The concept of spin-relaxation time is not, in general, applicable for non-degenerate pump-probe schemes, 
but the deviations from the exponential decay of the spin-polarization are confined to time scales of less than 100 fs, 
and can therefore only be resolved by a time resolution well below 100 fs. These results also demonstrate the limited 
validity of a single spin relaxation-time for holes. Strong p-type doping shortens relaxation times significantly due to 
the scattering from ionized donors. The calculated spin relaxation times support the assumption of very fast hole-spin 
relaxation in ferromagnetic GaMnAs at low temperatures. 
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